
DECUS NO. 


TITLE 


AUTHOR 


COMPANY 


DATE 


SOURCE LANGUAGE 


8-582 


RANDOM NUMBER GENERATOR ADAPTED FOR 
8K FORTRAN/SABR 


Adapted by: W. F. Haygood, Jr. 


3953 Starbrook Road 
Chesapeake, Virginia 

October 7, 1972 


8K SABR 


ATTENTION 


™l is a USER P r °9™- ° ther tha ? re 9 u ' r ' n 9 that it conform to submittal and review standards, 
no quality control has been imposed upon this program by DECUS. 

The DECUS Program Library is a clearing house only; it does not generate or test programs. No 
warranty, express or implied, is made by the contributor. Digital Equipment Computer Users 
Society or Digital Equipment Corporation as to the accuracy or functioning of the program or 
re ated material, and no responsibility is assumed by these parties in connection therewith. 















DECUS Program Library Write-up 


DECUS No. 8-582 


SUMMARY 

This is a Random Number Generator adapted from DECUS s-?s anH ^ 

£““ with DEC,S 8K ™™AN and 8K S»R pro^LSim 

identical lo DkJs’s-IJ? ° r P ° rU '’ n c< "‘ lne achie “ es •■«*«'»* 

The range -1 < x < +1 has been chosen as the output of the generator 
^‘thro'ug UlllT^Z ' DECUS + 5 - 25) outp ^s octal numbers in he range 

rrv-"™ 

the resui, can he ■ultiSi.S h/a sdiJahfl £25,*?**“ *" d ' Slred> 

first* I’oadPd'" initializcd *° the value 0000 0000 qm U when it ls 
first loaded. The user may initialize the generator to any desired value 
at any time by using a CALL to RAND I as explained below. Initializing 
the generator at the beginning of each Programming run will eliminate 

attained? ^ 9 nUmbCr ° f UmeS t0 advance to a Previously 

USAGE 

T ° 1 ?o^8K 1 FORTRANf nerat0r ^ Pr °* ram COntr °'’ USe the Allowing: 

CALL RANDI(J,K,L) 

m + K ! 3nd L are the decimal equivalents of the desired 36- 
bit initializing number (high* middle, and low order Parts 

random H° r example! if at . the end of a Program run. 

.. * lamed on page 2) contained the value 0200 1000 2000 fi » 

vafl F0RTR i N sta * cm '"< ““I" re-cs,abuses “ 

value at the beginning of the next Programming run: 

CALL RANDI(128, 512, 1024) 

for 8K SABR: 

CALL 3, RANDI 
ARG J 
ARG K 
ARG L 

where J, K, and L are the addresses of the high, middle, and low 
order parts of the desired 36-bit initializing number. For 

oaap P ?? ! lf + 9t t 5 e + ! nd ° f 3 program run ’ RANDOM (explained on 
P ?^ 2 H C °2 tain f d the value 0200 1000 20008, then the following 

L r r stabush ,his vaiue **• 

CALL 3, RANDI 
ARG (200 
ARG (1000 
ARG (2000 
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USAGE (continued) 


=SKSH^ 3 “-^»-iM= 


1) 

2 ) 

3) 

4) 

5) 


s^ d + ^V d S re \ S -° f * AND0M from the Linkinfi Loader Storage Map. 
+!* *5! Inst ™ction Field switches and the Switch Register to 
the address of RANDOM. Press LOAD ADDRESS. 

wo^! h %r“s“ C KPo“l!' r ,0 ‘ l ' Sired V3 ' UC ° f th€ hUh 0rter 

wo^d! h %re3“ C DEP0?jn er ,0 deSir ' d ValUe df * hc middlc »'- d ''' 

«o , rd. ,he pr“ss C SEPOS^! dr ‘° thC dCSlred ValUC ° f ,hc l0 ” order 


th€ n n c «3sa?rin3tn^! i meth0d * e '^ ina,eS the ne “ ssit * «♦ drodraomnd 
usinalhis mathod! S Cr ° firam - RANDI nced not bc called »" 


To call the generator, use the following call: 
for 8K FORTRAN: 

X=RANDU(I) 

where (I) can be any value. This value is ignored by the 
generator, so normally just (0) would be used thus: 

X=RANDU(0) 

Of course, the Random Number Generator call may be nested with 
other cal Is as in the following statement: 

X=S&RT(RANDU(0)) 

for 8K SABR: 

CALL 1, RANDU 
ARG 0 


To obtain random integers from the generato 
may be multiplied by a suitable constant, 
for 8K FORTRAN: 


r, the generator output 
For example: 


I=RANDU(0)*CONST 

where CONST is the constant. Note that the constant must be 

floating point since the generator output is less than 1 in 
magnitude. 


for 8K SABR: 

CALL 1, RANDU 
ARG 0 

CALL 1, FMP 
ARG CONST 


listing documentation consists of the SABR second pass 

listing, symbol table, and a sample FORTRAN Program showing the output 
of this generator for 320 random numbers. 
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//// 


//// 


//// 


NUMBER GENERATOR FOR 
FORTRAN/SABR SYSTEM 


//// 

/ 

/A RANDOM 
/DEC'S 8K 
/ 

/THE ACTUAL GENERATOR IS DECUS 5-25 
/ 

/ADAPTED FOR 8K FORTRAN/SABR BY 
/W. F. HAYGOOD JR., CWO-3, USCG 
/ 

/OCTOBER 7, 1972 
/ 

/THE ACTUAL 


CODING CONSISTS OF 2 PARTS: 

A SUBROUTINE WHICH INITIALIZES 
THE RANDOM NUMBER GENERATOR 
A FUNCTION WHICH IS THE 
GENERATOR ITSELF 


/ 1 ) 

/ 

/ 2 ) 

/ 

/ 

/THE GENERATOR MAY BE INITIALIZED TO ANY DESIRED 
/36-BIT VALUE BY MEANS OF THE FOLLOWING CALL: 

CALL RANDI (J,K,L) /FOR FORTRAN 

WHERE J, K, AND L ARE THE DECIMAL 
EQUIVALENTS OF THE HIGH, MIDDLE, 
AND LOW ORDER WORDS RESPECTIVELY. 
OR 

CALL 3* RANDI /FOR SABR 
ARG J 


/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/FOR 

/ 

/ 

/ 

/ 

//// 


ARG 

ARG 


K 

L 


WHERE J, K, AND L ARE THE ADDRESSES 
OF THE HIGH, MIDDLE, AND LOW ORDER 
PARTS OF THE DESIRED 36-BIT STARTING 
NUMBER. 

EXAMPLE, THE FORTRAN CALL: 

CALL RANDI(128, 512, 1024) 

WOULD SET THE GENERATOR TO THE OCTAL 
VALUE: 0200 1000 2000 (36 BITS USED) 


//// 


//// 


//// 
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//// 


//// 


//// 


//// 

/ 

/INITIALIZING THE GENERATOR WOULD ELIMINATE HAVING TO 
NU * BE * OF TIMES TO MAKE IT -22S5oJJ^ 1S21 ” 
/AFTER LOADING AND USING IT THE FIRST TIME# IT ALSO 

/?M L ?eItInG PROGR^S?” 0F " RAND ° M ” INPUT D * TA ™R USE 
/ 

/THE RANDOM NUMBER GENERATOR MAY BE CALLED AT ANY TTMF 
/WHETHER IT HAS BEEN INITIALIZED OR NOT. USE Jhe 

/following fortran statement to call it* 

/ X=RANDU<I> 

' WHERE I CAN BE AN* VALUE (USUALLY 0). 

/SINCE RANDU RETURNS NUMBERS IN THE RANGE -1 < X < +1 
/A SUITABLE CONSTANT MULTIPLIER CAN BE USED FOR 
/PRODUCING INTEGER NUMBERS THUS: 

/ J=RANDUC05*2048 

/ 

//// //// //// //// 
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//// 

//// //// 



/ 

ENTRY RANDI 
ENTRY RANDU 
ENTRY RANDOM 
DUMMY J 

DUMMY K 

DUMMY L 



random# 


0200 

0000 

HIGH# 

0 

0201 

0000 

MID# 

0 

0202 

0001 

LOW, 

J# 

K# 

1 

0203 

0000 

L# 

BLOCK 2 

0204 

0000 


0205 

0000 

RANDI, 

BLOCK 2 

0206 

0000 


0207 

4067 


TAD I RANDI 

0210 

0205 01 


0211 

1407 



0212 

2206 


INC RANDI# 

0213 

3203 


DC A J 

0214 

4067 


TAD I RANDI 

0215 

0205 01 


0216 

1407 



0217 

2206 


INC RANDI# 

0220 

3204 


DC A J# 

0221 

4067 


TAD I J 

0222 

0203 01 



0223 

1407 



0224 

3200 

/ 

OCA HIGH 



//// 

//// //// 


//// 

/generator initializer 

/GENERATOR ITSELF 


/RANDOM NUMBER FORMED HERE 


/DUMMY VARIABLES USED IN 
/INITIALIZATION PROCESS 


/GET ADDR OF J 


/GET VALUE OF J 


/HIGH ORDER INITIALIZED 
//// 
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0225 

4067 

0226 

0205 01 

0227 

1407 

0230 

2206 

0231 

3203 

0232 

4067 

0233 

0205 01 

0234 

1407 

0235 

2206 

0236 

3204 

0237 

4067 

0240 

0203 01 

0241 

1407 

0242 

3201 

0243 

4067 

0244 

0205 01 

0245 

1407 

0246 

2206 

0247 

3203 

0250 

4067 

0251 

0205 01 

0252 

1407 

0253 

2206 

0254 

3204 

0255 

4067 

0256 

0203 01 

0257 

1407 

0260 

3202 

0261 

4040 

0262 

0001 06 


/ 

//// 


//// //// //// 
TAD I RANDI 


INC RANDI# 
DCA K 

TAD I RANDI 


INC RANDI# 
DCA K# 

TAD I K 


DCA MID 
TAD I RANDI 


INC RANDI# 
DCA L 

TAD I RANDI 


INC RANDI# 
DCA L# 

TAD I L 


DCA LOW 
RETRN RANDI 


//// //// 


//// 

/GET ADDR OF K 


/GET VALUE OF K 

/MIDDLE ORDER INITIALIZED 
/GET ADDR OF L 


/GET VALUE OF L 
/LOW ORDER INITIALIZED 

//// 
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0263 

0000 

0264 

0000 

0265 

1202 

0266 

7004 

0267 

0377 

0270 

3361 

0271 

1201 

0272 

0376 

0273 

1361 

0274 

7006 

027 5 

7006 

0276 

7004 

0277 

3357 

0300 

1202 

0301 

0376 

0302 

7006 

0303 

7006 

0304 

7004 

0305 

1201 

0306 

3360 

0307 

7420 

0310 

2357 

0311 

7000 

0312 

1202 

0313 

7104 

0314 

7420 

0315 

2360 

0316 

7410 

0317 

2357 

0320 

7000 

0321 

7100 

0322 

1202 

0323 

3202 


/ 

//// 


//// //// //// 
/ 

RANDU, BLOCK 3 

TAD LOW 
RAL 

AND (7400 
DCA TEMP3 
TAD MID 
AND (177 
TAD TEMPS 
RTL 
RTL 
RAL 

DCA TEMP 
TAD LOW 
AND (177 
RTL 
RTL 
RAL 

TAD MID 
DCA TEMPI 
SNL 

INC TEMP 
NOP 

TAD LOW 
CLL RAL 
SNL 

INC TEMPI 
SKP 

INC TEMP 
NOP 
CLL 

TAD LOW 
DCA LOW 

//// //// 


//// 


/BEGIN DECUS 


//// 


-25 
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//// 

//// //// 



/ 


0324 

1201 


TAD MID 

0325 

7004 


RAL 

0326 

7420 


SNL 

0327 

2357 


INC TEMP 

0330 

7000 


NOP 

0331 

7100 


CLL 

0332 

1360 


TAD TEMPI 

0333 

3201 


DCA MID 

0334 

1200 


TAD HIGH 

0335 

7004 


RAL 

0336 

1200 


tad high 

0337 

1357 


TAD TEMP 

0340 

3200 


DC A HIGH 

0341 

2264 


INC RANDU# 

0342 

2264 


INC RANDU# 

0343 

1200 


TAD HIGH 

0344 

4033 


CALL 0, FLOT 

0345 

0004 06 



0346 

4033 


CALL 1, FDV 

0347 

0105 06 



0350 

6201 05 


ARG C 

0351 

0354 01 



0352 

4040 


RETRN RANDU 

0353 

0002 06 



0354 

2144 

C, 

2144 

0355 

0000 


0 

0356 

0000 


0 

0357 

0000 

TEMP, 

0 

0360 

0000 

TEMPI, 

0 

0361 

0000 

TEMP2, 

0 



/ 




//// 

//// //// 

0376 

0177 



0377 

7400 






END 


//// 


/END OF DECUS 5-25 
/SET RETURN POINTER 

/PUT RANDOM NUMBER IN FAC 
/PUT IN RANGE -1 TO +1 


/C»2048C10) 


//// 
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c 

0354 

FDV 

0000EXT 

FLOT 

0000EXT 

HIGH 

0200 

J 

0203 

K 

0203 

L 

0203 

LOW 

0202 

MID 

0201 

RANDI 

0205EXT 

RANDOM 

0200EXT 

RANDU 

0263EXT 

TEMP 

0357 

TEMPI 

0360 

TEMPS 

036 1 
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# 





c 


10 

20 

30 


SAMPLE RANDOM NUMBER GENERATOR 
GALL RANDI(18g, 512* 1024) 

DO 20 J*l, 40 
DO 10 1*1, 8 
X=»RANDU ( 0 ) 

WRITE (1,30) X, 

WRITE (1,30) 

STOP 

FORMAT (2H ,F7.4) 

END 


OUTPUT 


.1924 

.6060 

-.9062 

-.9404 

-.2344 

.8247 

.5122 

.8149 

.3809 

.1709 

.5303 

• 3691 

-.8481 

-.0156 

.0508 

.4307 

.7422 

.2549 

-.8901 

-.8916 

-.9375 

.7158 

.8867 

-.0610 

.4429 

.3569 

-.7695 

.4697 

.4214 

-.2075 

.7109 

.4111 

.4741 

-.5493 

.7305 

-.0303 

-.0659 

-.6694 

-.8042 

.8657 

-.3247 

.0542 

-.7134 

.6377 

.0674 

.7305 

-.5854 

.5215 

.9517 

.8828 

.1279 

• 1621 

-.8765 

.8989 

-.7104 

-.8530 

.0571 

-.7998 

-.5957 

.9912 

-.8555 

-.0386 

-.0869 

.0176 

-.5728 

.3101 

-.1685 

-.7275 

-.6860 

-.5293 

-.3828 

. 1294 

.9653 

-.0757 

-.7847 

.4243 

-.5342 

.9258 

.8135 

.7187 


-.1147 

-.1279 

.4961 

-.5459 

-.9590 

.8389 

.2622 

.2529 

-.4199 

-.0439 

-.5747 

-.7583 

-.4795 

-.7227 

. 1094 

.7935 

.8320 

.7119 

.6431 

-.1045 

.7148 

-.1392 

-.3657 

.3687 

.1382 

-.3701 

-.2739 

.1440 

.9453 

-.4482 

.0513 

.6196 

.4199 

-.5244 

-.7739 

-.3560 

.5591 

-.6069 

.4150 

.7109 

-.7700 

.9058 

.2305 

-.3428 

-.2417 

-.0093 

.3813 

-.3906 

.7153 

• 3618 

-.1963 

-.6226 

-.7358 

-.4912 

-.7422 

-.7607 

.6680 

-.7808 

-.7090 

.8384 

-.5483 

-.9292 

.8701 

-.9238 

.9976 

-.7915 

-.8652 

-.6299 

.9795 

.6543 

.2056 

.0825 

.8403 

-.2632 

-.4097 

-.2632 

.3447 

-.2505 

.9736 

- .6128 


.2461 

.6421 

.2422 

.3799 

-.3535 

.3428 

-.8613 

.5693 

-.5015 

-.6011 

.6050 

.4683 

-.0381 

.2886 

-.2397 

-.5654 

.7651 

.1982 

-.4331 

. 3 540 

.7129 

-.4580 

-.5146 

-.3911 

.5195 

.4619 

-.6880 

.5894 

.7847 

.7573 

-.7612 

-.1309 

-.9434 

-.9268 

• 8120 

.0898 

-.6909 

-.6689 

.5127 

.6914 

.3472 

-.0547 

-.1489 

.2061 

.1021 

.7095 

.2075 

.7759 

-.2842 

-.9492 

• 0122 

-.3105 

-.3437 

.3716 

.0977 

-.5537 

-.7129 

. 7 646 

—.6050 

.8384 

-.6572 

.4316 

.6084 

-.0215 

.2554 

.6694 

-.0103 

-.3774 

-.9077 

.6802 

.6274 

-.9644 

.8403 

-.5757 

.0903 

.9243 

-.6216 

.5391 

-.4580 

.7793 


-.3813 

-.0537 

.0820 

-.9131 

-.7793 

.2524 

-.8496 

-.2080 

.8892 

.7583 

-.6528 

-.1177 

.0581 

-.2349 

.7461 

-.4209 

.2871 

-.0483 

.0039 

.8525 

.8179 

-.9561 

.2686 

-.8540 

.0781 

.3252 

-.2896 

.9722 

-.5352 

-.0146 

• 0493 

-.5103 

.9116 

-.1748 

-.7881 

.4766 

.1855 

-.8530 

-.4834 

.8901 

.5288 

-.3228 

• 5601 

-.4790 

-.6792 

-.4468 

.7695 

-.3511 

.8447 

-.3755 

.0107 

.8730 

-.6929 

.5117 

-.2187 

-.3159 

.9863 

-.9502 

.4575 

-.7866 

.4873 

-.9473 

.3779 

.4741 

-.2993 

.1924 

-.1904 

.2695 

.2324 

-.7124 

.5488 

-.0132 

.9653 

.9368 

.7153 

-.0132 

.8096 

.0195 

.7803 

-.3188 
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